Letters Author Response: Passing–Bablok Regression Is Inappropriate for Assessing Association Between Structure and Function

We thank Marin-Franch for his comments on our recent paper. If some aspects of the statistical analysis were not as clear as they perhaps could have been, we are grateful for the opportunity to address them here. The first concern was that our conclusion was in conflict with previous findings. Hood et al. assumed that there is no relationship between retinal nerve fiber layer thickness and visual field sensitivity in healthy observers; this was not a finding, but a working assumption. This point was made in our paper. On the other hand, Wollstein et al. suggested that substantial structural loss must occur in the retina before functional damage can be detected using current testing methods; however, other published data support the notion that retinal structure and visual function are associated in healthy eyes. Marin-Franch’s second, and principal, concern is our use of Passing–Bablok regression analysis to test the null hypothesis that there was no association between retinal structure and visual function; he instead suggests that more insight would be gained from applying ordinary least squares (OLS) linear regression or Spearman’s correlation. Firstly, we would like to point out that Marin-Franch has mistakenly applied the ‘‘method-equality’’ analysis outlined in Passing and Bablok in his simulations, rather than the transformation method, as outlined by Bablok et al., which we used in our analysis. Critically, for the derivation of the transformation regression procedure in Bablok et al., the authors state that it is necessary to start with the assumption of a positive correlation, as determined by Kendall’s correlation, which is a similar, and an alternative, rank-based correlation method to Spearman’s correlation, as proposed by Marin-Franch. When we performed Kendall’s correlation test on our combined data, we did indeed find a significant association. Thus we felt it appropriate to use this technique. We regret that this point was not made in our paper, but it is clearly stated in the correct reference. In any case, inspection of Figures 2 through 4 in our paper reveals a striking similarity in the slopes of Passing–Bablok regression lines between hemifields and between healthy and glaucoma data, for most comparisons. Is it better to use OLS here? We would argue that it is not. Warton et al. point out that it is inappropriate to use OLS when the variance in the xand y-variables differs and when noise is not Gaussian. Furthermore, they also state that an errors-in-variables method like major axis regression or standardized major axis regression should be used if one wishes to determine the slope parameter after establishing a significant association, which, as already pointed out, is an assumption of the method described by Bablok et al. Furthermore, part of our motivation for performing Passing– Bablok regression on our data came from analysis of 200,000 simulated data sets with Deming regression and Passing– Bablok regression. These simulated data sets comprised total deviation values (standard automated perimetry) and retinal nerve fiber layer thickness (optical coherence tomography), based on the data of Hood et al., and the aim was to determine how well each regression technique performed at estimating the true underlying slope in the simulated data. Marin-Franch and colleagues previously showed, also using simulations, that estimating linear relationships between structure and function measurements in glaucoma may be improved by the use of Deming regression compared with OLS, major-axis, and standardized major-axis regression techniques. However, unlike Passing–Bablok regression, Deming regression requires the investigator to input the ratio of variance in the two measurements of interest from population data. In practice, this information is often unavailable, as was the case for peripheral grating resolution acuity in our study; moreover, the variance ratio may change